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Abstract 

Ginzburg-Landau functional is derived for a system possessing both chiral and diquark 
condensates. Anderson theorem for such a system is formulated and proved. 

1 Introduction 

During the last couple of years color superconductivity became a compelling topic in QCD - see e.g. 
the review papers ||^, ^ . The present work is aimed to shed some new light on the region of mixed 
condensation in which both chiral and color symmetries are broken. For definiteness we consider 
the case of two flavors u and d, the corresponding color superconducting phase is called 2SC |jl|, 
The point at issue is the interplay of the chiral condensate ip and the diquark condensate A in 
the interval of the chemical potential fi where they possibly coexist. Necessary to mention that 
as model calculations show Q-Q the very presence of such a region is questionable. Instead of 
coexistence one rather observes a strong competition between ip and A in the sense that where one 
condensate is nonzero the other vanishes. This conclusion has been drawn as a result of numerical 
analysis of the complicated set of two coupled gap equations for ip and A ||, |^ (other authors take 
it as a plausible apriori assumption Q ) . The question has been also addressed within the random 
matrix approach |j] with the result that for single-gluon exchange the coexistence region of the 
two condensates is absent. Evidently such a situation calls for simplification. Arguments presented 
below are aimed to show that as soon as (with /i increasing) the diquark condensate A develops the 
role of the chiral condensate ip is prescribed by what may be called the Anderson theorem ^ . In 
BCS theory this theorem states that nonmagnetic impurity do not influence the thermodynamic 
properties of a superconductor (in particular the values of the gap and of the critical temperature) 
in linear approximation in impurity concentration. In case of color superconductor the role of 
impurity is taken by chiral condensate. This means that the system of the gap equations for ip 
and A actually decouples and with high accuracy one can calculate A neglecting the influence of 
ip. One general remark is probably needed. Namely, the content of the "true" Anderson theorem 
in BCS is physically deeper than the statements formulated below. This point has to be discussed 
within comparative study of BCS theory and color superconductivity which is beyond our present 
goals (see [1,2]). 

In Section 2 we present and discuss the general expression for thermodynamic potential which 
has the structure common to three approaches to color superconductivity, namely NJL, one gluon 
exchange, and instanton ensemble. In Section 3 we derive the Ginzburg-Landau functional for 
chiral condensate, while in Section 4 the same is done for that part of the potential which depends 
on both condensates - the diquark and the chiral ones. The main result of this Section is that the 
contribution of the chiral condensate to the " mixed" part of thermodynamic potential may appear 
only in orders higher than second one. In Section 5 we show that the same is also true at zero 
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temperature. Section 5 is devoted to the formulation and discussion of what we call the Anderson 
theorem for color superconductor. 



2 Partition Function and Thermodynamic Potential 

We start with the QCD Euclidean partition function 

Z ^ I DAD^pDtPexpi-S), (1) 



where 

^=\J = j ^Pi-tj^D^ - zm + *M74)V'rf^x. (2) 

In (H) color and flavour indices are suppressed, Nf = 2, Nc = 3, and the chemical potential 
fi is introduced. Performing integration over the gauge fields one gets effective fermion action in 
terms of cluster expansion 

Z = J D^PDiPexpi- J d'^xLn- Seff), (3) 

with Lq — ip{~i'j^d^ — im + iij,j4)ip and effective action Seff — J2'^=2 ;7r(i£"))' where 6 ~ 
J d'^x'ijj{x)g'yfj,A'^{x)t°'ip{x) and double brackets denote irreducible cumulants ||lO[. Since we con- 
sider only u and d quarks the current quark mass m will be neglected. In future we plan to give 
up this approximation especially in view of the numerical results obtained in which show that 
there is some influence of m on the phase structure near the transition point. One should also 
keep in mind that integration over gauge fields when passing from (|l|) to (^) is in no way a trivial 
operation. It certainly deserves dedicated investigation which is outside the scope of the present 
work (for some guidlines see ijlTI). 

To proceed further some simplification of (||) is needed. Being quite general, (||) gives rise to 
several commonly used models. First step is to keep in Seff only the lowest four-quark interaction. 
Routine manipulations then lead to the one gluon exchange model JT2| which in turn displays 
the bulk properties of the color superconductivity phenomena [|]. One gluon exchange yields 
surprisingly interesting results at very high density ^ . Less trivial but very transparent approach 
psf enables to recast (||) into the delute instanton gas model which is successfully used in color 
surerconductivity studies ||, 1). We note that derivation of the instanton model presented in 
pagoes beyond zero modes approximation. The importance of higher modes was demonstrated 
in ]15|, [l^ . Simplifying things even further [p3| one arrives at the N JL models which as also quite 
suitable to describe color superconductivity ^ . 

A lesson from comprehensive studies of color superconductivity within different models listed 
above is that the main results are very similar ^ . The key technical points leading from chosen 
Lagrangian to the effective action are similar as well. Use is made of the bosonization of the fields 
{'tjjip) and {ipip) Hi- As a result the following general expression in the Nambu-Gorkov basis of 
eight component fields emerges: 

Seff= [d^x{^ + ^-'-trln(^%, T ■ ,T^. ff'^+^^^-^^^^^ )\. (4) 

Let us present necessary explanations to (^). We use the following representation 
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C = 7274, C-'i^.C^-^l,dl^t^ = -^,. (6) 

Effective action (||) is written in terms of the two condensates - the chiral one ip and the super- 
conducting diquark one A. Operators $ and A in determine the matrix structure of these 
condensates 

Af/ = ea/33e»jC75A = $A, <^f/ = S^pS^^e = Kip, (7) 

where indices a,/3 correspond to color, and i,j - to flavor. The matrix structure defined by (|^) is 
the simplest possible. For example, if one starts with the single gluon exchange model and performs 
Fierz transformation with respect to Lorentz, color and flavor indices, one arrives to several other 
diquark condensates in addition to (|^), e.g. the negative parity one with <&' — Sa/33£ijC, etc., see 
[0, The same is true for the matrix A. The minimal "improvement" of A in line with chiral 
invariance requirement would be A' = Ta^o-s + *75''"a=o-3 while in only the a = part of the 
first term in kept. In other words from the two fields {a, tt) only a is kept. The two coupling 
constants gf and have the dimension of m~^. In the one gluon exchange inspired model they are 
related to each other but in most studies they are considered as independent |l[ . Regularization of 
the integrals involving the chiral gap is performed either by smooth formfactor or by cutoff. The 
sensitivity to the regularization procedure is unessential and to avoid complicated equations we 
do not introduce this regularization explicitly. Integrals involving the superconducting gap A are 
regularized at the Debye frequency -see below. 

Factorizing from (j^ the four-volume and making use of the identity 

trln(^^^T '^+^ =tr\nAA^ + trln{l + A-^R.{A'^)-^R+} (8) 

one arrives at the following expression for thermodynamic potential Q{ip, A; jj,, T) 

A+A 
4ff! ^ 



fl{p, A : fj,,T) — ^j— 2 + ——2 tr ln(-p + iipA - i^i^i)- 



1 , AA+$<i>+ 1 

- -tr\n{l + — — — — ttL (9) 

where p — Pfj.jfj,,p± — {pk,P4 ^ip)- This form on was first derived in j^. 

The remaining part of the paper will be mostly devoted to the investigation of in the 
Ginzburg-Landau region, i.e. to the expansion of 17 in powers of A and p at temperatures close 
to the critical one. Obviously the system under consideration possesses two different critical tem- 
peratures - the chiral temperature and the diquark one Ta. The situation simplifies if one 
focuses on the transition region in the (T, p) plane where the two condensates and A compete. 
Numerical calculations show |^ that in the corresponding interval of /i the slope ^ of the phase 
transition line is small and a notion of a common critical temperature Tc makes sense. At T > Tc 
both condensates ip and A vanish. Numerically Tc is about 30-40 MeV while the corresponding 
values of /i lie in the interval of 300-400 MeV . This is the region of the phase diagram which 
we will be interested in. Consideration will be restricted to the case of homogeneous condensates 
p and A. Thus the gradient terms will be absent. 

We start the derivation with the expansion of the third term in in powers of p). 

3 Chiral Condensate Transition Line 

Consider the "tp-part" of the partition function (^), i.e. the functional 

2 

n^{ip;fi,T) = — ^ -irln(-p + z(^A-z^74). (10) 
451 
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The equation of Q^p proceeds along the standard hnes One has 

d^p 1 



)3 [(2n + 1)ttT + +p^ + if' 



2 ■ 



(11) 



where Nc — 3, Nf — 2, and the summation is over the fermionic Matsubara modes. After perform- 
ing the sum and dip'^ integration the resuh for fl^p reads 

= ^ - j dpp^{ln2che^-^ +In2ch9''+'^), (12) 



Next we expand Q^p in powers of ip in the vicinity of T^. Expansion of the thermodynamic 
potential containing the second and the fourth order terms is called Ginzburg-Landau functional. 
The first order term is absent because of the mean- field stationarity condition dflip/d(p — 0. The 
discussion of the interesting physics associated with the tricritical point and sixth order term is 
out of the scope of the present paper, see ||l|, ^, |l^, |2^. 

With the expression ( [T^ for the thermodynamic potential at hands the calculation of the 
derivatives d^ft^/dip^ and d'^ilip/dip'^ is rather straightforward. One has 



2A(T) - ^ 



^ _1 NrNf 



dip 

Denoting 6'^^^ at T = Tc as oi"^^ one writes 



2.9? 27r2 



dppithO'--^ +the'-+^). (14) 



Substitution of ( |l5| ) into (|4|) leads after simple integration to the result 

^(r, = iw,"(^)-T,"(I^). (16) 

In arriving to (|l^ use was also made of the gap equation dVl^/dip = which yields cancel- 
lation of the 1/2^2 term. Strictly speaking the gap equation is dft/dip — 0, where SI is the full 
thermodynamic potential (^ but as we shall see the corresponding corrections are small. Note 
that A < at T < Tc in line with the general theory of phase transitions. 

Calculation of the fourth order derivative is almost as simple. As a preliminary remark we 
remind that for Nf ~ 2 chiral symmetry restoration at fi = proceeds via second order phase 
transition. On the other hand at T = but with fi increasing the restoration of chiral symmetry 
occurs via first order phase transition ^, ^ ^ . Correspondingly in the first case the coefficient 
of the quartic term has to be positive while in the second case -negative. The sign is changed at 
the tricritical point ||l|, ^, |l^, ^ at which the sixth order term is important. Consideration of this 
point is beyond the scope of the present paper. Another remark is that in the NJL model the order 
of the phase transition depends upon the value of the cutoff |2^. We shall immediately see that 
indeed the sign of the fourth order term depends on the value of the chemical-potential /x and of 
the cutoff value A. 

The nonvanishing at ip = part of the fourth order derivative reads 

, d-^n^ 3NcNf f , r th0^-'> + the^+^ i / i i M 
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Simple manipulations yield 



As expected the sign of B depends on the values of the chemical potential ^ and cutoff A. Having 
in mind to demonstrate explicity the dependence on the cutoff we may resort to a rough estimate, 
namely 

A A /"^ , Ini 
■ In — - / dt- 








/■^/^^ _ tht 
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TT 



Tc Tq ,/n ch^t 
t/i-f In-l+ln^^lnl^, (19) 

where 7 = e^, C = 0.5777... Typical cutoff value is A - 800 McV |, |, |2l). 

Collecting pieces together we write down the Ginzburg-Landau functional for the chiral part of 
the thermodynamic potential {Nc = 3, TV/ = 2); 

(20) 

We note that the above treatment was oversimplified in several points, in particular we assumed 
that the coupling gf was temperature independent, the quark current masses were set to zero, etc. 
These omissions important by themselves are hardly crucial for our main purpose which is the 
formulation and proof of the Anderson theorem. 




4 Ginzburg-Landau potential for diquarks 

Now we turn to the last term in Eq. (|9|). Expanding the logarithm in powers of we directly 
arrive at the desired Ginzburg-Landau functional. The term proportional to reads 

45i 2'"(p+-z^A)(p_+z^A)- ^''^ 



Inverting the denominator and performing the trace over discrete indices one gets {Nc = 3) 

o(2) - SA^tV f '^'P (p1 + p^ + m^ + ^^) 

where 

R^pt + 2p2(p2 +^^+ (^2) ^ (p2 _ ^2 ^ ^2)2^ (23) 

and the sum (|2^) is taken over fermionic Matsubara modes corresponding to p^. Next step dif- 
fers from the calculation of the (p'^ term presented in the previous section. Namely we assume 
that considering the Landau-Ginzburg region for the diquark condensate A one can neglect the 
contribution of antiparticles. This means that in the representation 



R 




(24) 
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only the first term is kept. Then summation in (|2^) yields 

.2 



A2 ri^n 1 „ 

where 

£.=p-i-t, q = Vl^T^ - fJ; (26) 

and lud is the Debye frequency while A^(0) is the density of states at the Fermi surface 



2_ / ^dp\ ^ V 

[P di) p TT^ 



N{0)^^[p'-fA (27) 



The density of states ( p7| ) is four times larger than the BCS theory factor /i^/27r^. This is due to 
color and flavor degrees of freedom, in 2SC there are red and green u and d quarks. 
Next we expand the integrand in ( p5|) in powers of (T — Tc) and ^p. This gives 

with y = ^/2T'. The last term drops upon integration over the interval [— a;j)/2Tc, cjj)/2Te]. 
Integration results also in cancellation of the first term with the term A^/4(7| in ( pl| ) (because of 
the gap equation or equivalently from the definition of the critical temperature). Thus 

= A'N{Q) [^^) ■ (29) 

We see that up to cp'^ there is no dependence on the chiral condensate, therefore we dropped the 
index ip in (|2^). At this point one may notice the analogy with the Anderson theorem in BCS 
theory. 

Next we consider the term proportional to in the expansion of the logarithm in (|^). In 
this term the ip dependence will be neglected from the beginning. There is no need to present the 
details of the calculations which are simple and essentially repeat the derivation of the A^ term 
with only two new technical points which are 

I d ^ 1 f°°dyfthy\'_ 7C(3) 



The result for fi^^'* reads 



(31) 

Thus the development of ^l^if in powers of A and p> reads 

f]A = -A2iV(0)(l- — -^^^A^l. (32) 

We note that with realistic values of A and /i the numerical value of Qaip given by ( |3^ is 
in perfect agreement with our estimate presented in . 

Now we are tooled to formulate conclusions concerning the interplay of the two condensates. 
However we first turn from the Ginsburg-Landau region to the T = case in order to demonstrate 
that the version of the Anderson theorem also holds at T = 0. 
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5 Renormalized Potential at T = 



At T = the formal development of the thermodynamic potential (|^) in powers of A and is not 
appropriate. Instead we calculate the potential explicitly. The "(/s-part" of it was evaluated in Sec. 
3 exactly so we may just set T = in (^2|) and get 



= 0) = |i - ^ / dvv\^^^Tv^^-{ii - V5MV)^?(/^ - V5^T^)}- (33) 



92 

Next we turn to the last term in 

f^A>p = — ^ - -ir In n + — -— — \ . (34 



Here use can be made of the standard trick; first one takes A0a(^/9A, then performs integration 
over dp4 which corresponds to T = 0, and then reconstructs f2Ai^(T = 0) by integration back over 
A. The result reads 

f^A^(r = 0) = ^-^ 1 dp/{\/(v^^ + ^-A*)2 + A2 + y(^^^T^-^)2 + A^}. (35) 

As before the second term corresponding to antiparticles will be neglected. In order to avoid 
the onset of " irrelevant" terms of higher order in A we renormalize f^Aip (2^ = 0) at certain scale 
fj, pl] . The renormalized coupling constant is defined as 



5A2 



A = K. (36) 



The renormalized potential reads 

i2 r A 2 A 2, ^2 



A2 r A2 A 2 2 



(37) 



where 



i?a = {(vV+^-/^)' a = A,K. (38) 
Performing integration and expansion in ip one gets 

Ok. = ^-iV(0)A2g-ln^)~A^(0)|;. (39) 

We see that correction of the order ip'^ is absent. 



6 Assembling the Pieces 

The main result of the paper is the Ginzburg-Landau functional with both condensates included. 
It reads 

9 4 9 f /"^/'^'^ , tht f cht \ ] 



\ Tc J 2Tr^ ] Jq i I cht — ch^- I j 
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The first two terms stem from tlie "cliiral" part of the thermodynamic potential while the 
"mixed" part resulted in the last two terms with contribution suppressed (only terms of the 
order higher than tp^ may appear). Comparing the coefficients of the first two terms with those of 
the second two one notices that the contribution of the chiral condensate is damped by a factor 
T2//x2 <c 1{Tc ~ 40 MeV, ^J^ ~ 400 MeV). 

Therefore we arrive to the following two conclusions; 

(i) The system of gap equations for chiral and diquark condensates actually decouples. 

(ii) As soon as diquark condensate is formed the contribution of the chiral condensate to 
thermodynamic quantities becomes strongly suppresses. 

In analogy with the BCS theory these statements form what may be somewhat loosely called 
the Anderson theorem. 

Its validity is supported not only by numerical calculations within several modes |3|-|^ but also 
within the random matrix approach based solely on the symmetry of the system. 

Several assumptions and simplifications made on the way to the above conclusions were outlined 
in the text. They are the following. The current quark masses were neglected and the same critical 
temperature was taken for chiral and diquark condensates. The first assumption is a technical one 
and was made in order to avoid lengthy equations. It will be lifted in the next publication. The 
equality of the critical temperatures and Ta in the vicinity of the ^ value where the diquark 
condensate arises was demonstrated by numerical calculations in [||. This assumption is neither 
crucial and our final result (^) is easily generalized to the case of different T^p and Ta. The 
expression for the thermodynamic potential is quite general and is valid for a variety of the 
interaction models and different values of the coupling constants. The central step is the expansion 
of the logarithm in (^) up to and terms. This procedure is valid near the transition point 
but beyond the region of strong fluctuations. 

The suppression of the chiral condensate by the factor T^j [P' in the final expression (4^) for 
is due to the /x^ factor in the density of quark states at finite chemical potential (see (p7)). 

Finally we want to stress the point which we consider most important for future amendment. 
This is the investigation of the role played by gluon fields Nonperturbative gluon fields 

introduce into the subject a new quantity - the correlation length Tg. It might be even speculated 
that the term A'^ in the Ginzburg-Landau functional will appear. 

The author is grateful to Yu. A.Simonov, V.I.Shevchenko and N.O.Agasian for discussions. 
Support from grants RFFI-00-02- 17836, RFFI- 00-15-96786 and INTAS CALL 2000 project # 110 
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